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Abstract 

A solution of the Einstein's equations that represents the superposition of a 
Schwarszchild black hole with both quadrupolar and octopolar terms describ- 
ing a halo is exhibited. We show that this solution, in the Newtonian limit, 
is an analog to the well known Henon-Heiles potential. The integrability of 
orbits of test particles moving around a black hole representing the galactic 
center is studied and bounded zones of chaotic behavior are found. 
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There are two main lines of research of chaotic behavior in General Relativity: one deals 
with chaos associated to the existence of generic singularities in homogeneous cosmological 
models, as in Bianchi IX model In this case, Einstein's equations themselves - although 
suitably reduced to a simple toy model - form the dynamical system being considered. 
Invariance of the relativistic theory against arbitrary coordinate changes, particularly time 
changes, challenges the proper concept of chaos in this approach . The other line assumes a 
given geometry and looks for chaotic behavior of geodesic motion in this background. In this 
case, the geometry can be taken as either an approximate or an exact solution of Einstein's 
equations. Examples of chaotic geodesic motion in exact geometries are considered in . 
On the other hand, a linearized treatment of the geometry is made in and [0], to model 
the interaction of a black hole with a weak gravitational wave and a weak gravitational 
dipole moment, respectively. 

Firstly in this letter we consider - with some astronomical motivation - an exact solution 
of Einstein's equations that represents a Schwarszchild black hole together with quadrupolar 
and octopolar contributions of arbitrary strength. Since the multipolar contribution can be 
switched off, this solution is considered to represent a perturbation of a black hole. Secondly, 
we address the question whether this exact perturbation breaks the integrability of test 
motions around the black hole. Some previous results [^,0 deal with linearized perturbations 
of the Schwarszchild background which, although interesting, make the perturbation "strong 
enough" to push it out of the full theory. This is an important limitation since we know 
that chaos is rather the rule than the exception in the general context of dynamical systems 
and that only a small change in the dynamics could drastically alter the motion portrait. 

These facts motivate us to look for a solution of vacuum Einstein's equations whose 
additional controlable parameters are linked to multipolar structures, originated in the 
model from exterior matter distributions. To this end, our search has been inspired by 
the well known Henon-Heiles potential (HHP) - classical in astronomy and also in dynami- 



cal systems- one of the most studied paradigms of chaotic behavior PHin]. It was originally 



intended to model a Newtonian axially symmetric galactic potential. In fact, we have found 



a solution that is a general relativistic analog of HHP and moreover seems to be a less ide- 
alized model for an axially symmetric galaxy with a massive core plus an exterior halo of 
quadrupoles and octopoles. 

Finally, we shall see that the quadrupolar terms only are not sufficient to cause chaotic 
geodesic motion. The chaotic behavior occurs only due to the presence of the octopolar 
terms, which break the reflection symmetry of the potential about the middle plane of the 
galaxy, in the same manner that the cubic terms in the HHP give rise to the nonintegrability 
of the motion. 

We start with the Weyl metric describing static axially symmetric space-times 

ds' = e^'c^dt^ - e-'-'ie^^dz^ + dr^) + r^dcj^^ (1) 

V and 7 are functions of r and z only. The Einstein's equations reduce in this case to the usual 
Laplace equation v^rr+T-'^r/'^+T-'^zz = and the quadrature d'~^ = r[{h'rY — {h'z)^]dr+2rh'^r^zdz. 
Despite the simplicity of the Einstein's equations in the coordinates r and z, they are rather 
deceiving better coordinates are the prolate spheroidal u and v defined as 

u = ^[^r^ + {z + Lf + + {z - Lf], m > 1, 

V = ^[^r^ + {z + Lf - Jr^ + {z-L)% -l<v<l, (2) 

where L is an arbitrary constant. In these coordinates Laplace's equation can be solved in 
terms of standard Legendre polynomials. In particular, we are interested in the solution 

p{u,v) = aoQo{u)+b2P2{u)P2{v)+bsP3{u)Ps{v). (3) 

The first term is the monopole one, which is related to the Schwarszchild solution with 
Oo = — 1 and L = GM/c^ = Rq. The remaining terms represent the multipolar structure of 
a halo. We impose the additional condition of elementary fiatness over r = for \z\ > L, 
i.e. , 7(m,w) = in this region of the z axis, thus eliminating conical singularities there. 

The integration of 7 is guaranteed by the fact that u is a solution of Laplace's equation. 
We find 
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u — 1 

2i/ = (1 + e) log(-^) + z/q(m, v) + uo{u, v), 
-1 

27 = (1 + ef \og{— -) + 7q(m, v) + 7o(m, v) + 7qo(m, v), (4) 



(5) 



where e is an arbitrary constant and 
UQ = {Q/3){3u'-l){3v'-l), 
uo = iO/5)uv{5u^ - 3){5v^ - 3), 

7q = -4Q(1 + e)M(l - v^) + (QV2)[9mV - lOu^v^ - 10m V 

+ 12M2t,2 + ^4 ^ ^4 _ _ 2^2 ^ 

70 = -20(1 + e)[3u^v - 3mV + t;^ - ft; + |] + 202[f - IIImV 

_ 117^4^6 + |^6^2 ^ |^2^6 + 1|9^4^4 _ |Z^4^2 _ |Z^2^4 ^ |1^2^2 

_ 3 6 _ 3 6 , 27 4 , 27 4 _ ^ 2 _ ^ 2 , ^1 
8" 8 ^ 40 ^ 40 200 200 ^ 200 J' 

7go = 2QO[9u^v^ - 12u^v^ + 3u^v - 12u^v^ + fu^v^ - fuH + 3uv^ - fuv^ + ^uv]. 
Q and O are respectively the quadrupole and the octopole strengths. The terms proportional 
to Q and O in 7q and 70, respectively, represent nonlinear interactions between the black 
hole and the multipoles. For e = and L = Rqws have the Schwarszchild solution perturbed 
by quadrupoles and octopoles and, if e = — 1, we are left only with the multipolar structure. 
Note that only the odd terms in v break the reflection symmetry about the plane z = 
and all them are linked only to the octopole strength. The particular case e = O = and 
L = i?o is considered in ||TI|| . 

To see that the multipolar structure in the solution (|) and @ is a general relativistic 
analog of HHP, we write the solution in the usual Schwarszchild coordinates (t, p, 9, 0) 



z = {p — L) cos 9, r = p{p — 2L) sm9. (6) 
For e = and L = i?o in @) plus we get 

ds^ = (1 - ^)e(-e+-o)c2rft2 _ e{7Q+7o+7Qo-.Q-.o)r(i _ '^y^dp^ + p^d9^] 
P P 

+ 6-'^"^+"°^ p^ sin^ 9d(f)\ (7) 

In the Newtonian approximation only gu in (|^ is important. We expand it to the first 
order in the strengths Q and O and identify a Newtonian potential $Ar through 



gtt = l + ^$7v • (8) 
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Remember that (|]) are valid only for p > 2Rq. If we suppose that there exists some region 
between the horizon event p = 2Rq of the galactic core and the external multipolar structure 
simulating the halo, in which the Schwarszchild coordinates are approximately the Euclidean 
spherical ones, then there we have z ^ p cos 6 and p"^ ^ z"^ + R^ with R being the radial 
cylindrical coordinate. Expanding and using this approximation, we arrive at 

^N = -^ + 4^A^z'-R' + fip,z)] + '^[2z'-3zR' + 9ip,z)] , (9) 



p 2Rt '''' " 2Rf 



where 



f(34-l)(12f -14 + 4^; 
3 -Ro P 



/(p,z) = -^(3--l)(12^-14 + 4 

5 "'"p^ '^'^"'^ Rq " ■ — ^ 



^(p, z) = -^z{5-^ - 3)(25^ - 42 + 26^ - 4^). (10) 



For comparison we quote HHP from ^ 

^BH = l[z' + R' + 2zR'-lz% (11) 

The quadrupolar and octopolar terms explicitly shown in (|^) are qualitatively the same 
as those appearing in HHP, with one difference that will be of importance in the latter 
discussion: in HHP the quadrupolar terms constitute an (isotropic) bidimensional oscilator, 
while in our case we have a saddle point. Once HHP is considered to be by far more 
interesting by its dynamical richness allied to simplicity than by actually modeling any 
concrete astronomical situation, the similarity above justifies the parallelism. Moreover, the 
quadratic and cubic terms in (y) are indubitably in the extended Henon-Heiles family of 
potentials 0. 

In Fig. |1| we show Poincare sections across the plane ^ = in Weyl coordinates for 
some time-like geodesic orbits for the metric (^ and (^. In Fig. 1(a) we show a section 
for the unperturbed Schwarszchild background (e = 0, Q = O = 0), where the complete 
integrability is seen. In Fig. 1(b) we switch on the octopolar terms over Fig. 1(a) (e = 0, 



Q = and O ^ 0). This breaks the reflection symmetry about the middle plane and hence 
also the integrability of the motion. We also have checked that this chaotic manifestation is 
not explained by the signal structure of the eigenvalues of the curvature tensor as proposed 
in [Q. We see that the octopole terms are a powerful source of chaos since the strength we 
have used is only O = 1.5 x 10^^ in nondimensional units. In Fig. 1(c) we switch on the 
quadrupolar terms over Fig. 1(a) (e = 0, Q 7^ and O = 0). We observe that integrability is 
preserved, yet in a more involved toroidal structure in phase space. Finally, Fig. 1(d) exhibits 
chaotic behavior due to the full perturbation of the Schwarszchild background (e = 0, Q 7^ 
and 0^0). We note that although the quadrupolar terms do not cause chaos, their presence 
only with strength Q = 2.2 x 10"^ does the system more sensible to the octopolar terms. 

By switching the black hole off (e = — 1) - with or without the octopoles - we were not 
able to flnd regions of bounded motion that include points outside the middle plane. This 
lack of nonplanar bounded regions when solely the multipolar halo is taken into account is 
probably related to the repulsive quadractic 2;-term present in the expansion (|^). 
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FIGURES 

FIG. 1. Poincare sections (in Weyl coordinates) in the plane 2; = for timelike geodesic orbits 
in different geometries. All test particles (with mass /i) have energy E = and angular 

momentum L = 3.8RofJ-c. In all figures e = and the nondimensional variables are defined by 
X = r/Ro and dT = ds/Ro'. (a) unperturbed Schwarszchild metric {Q = O = 0). (b) Schwarszchild 
metric perturbed by octopolar terms (Q = 0,0 = 1.5 x 10^''). (c) Schwarszchild metric perturbed 
by quadrupolar terms (Q = 2.2 x 10~^, = 0). (d) Schwarszchild metric plus the full perturbation 
(Q = 2.2 X 10-^,0 = 1.5 X 10-'^). 
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